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ABSTRACT 

We present a scenario for the formation of disks which explains not only the 
properties of normal galaxies, but the properties of the population of low surface 
brightness galaxies (LSBs) as well. We use a gravitationally self-consistent 
model for disk collapse to calculate the observable properties of disk galaxies 
as a function of mass and angular momentum of the initial protogalaxy. The 
model naturally produces both smooth, asymptotically flat rotation curves and 
exponential surface brightness profiles over many disk scale lengths. 

In this scenario, low mass and/or high angular momentum halos naturally 
form low baryonic surface density disks, which will tend to be low surface 
brightness. Theoretical and numerical calculations suggest galaxy halos should 
form with a wide range of mass and angular momenta, and thus, the disks which 
form within these halos should have a wide range of surface brightnesses and 
scale lengths. We use theoretical predictions for the distribution of halo masses 
and angular momenta to explicitly calculate the expected number density of disk 
galaxies as a function of central surface brightness and disk scale length. The 
resulting distribution is compared to the observed properties of galactic disks, 
and is shown to explain the range of observed disk properties, including the 
cutoff in the maximum disk scale length as a function of surface brightness. We 
also show that disk instabilities explain the observed lack of high surface density 
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disks. The calculated distribution of disk properties also suggests that there 
are large numbers of galaxies which remain undetected due to biases against 
galaxies with either low surface brightness or small scale length. We quantify 
this by calculating the difference between the intrinsic luminosity function and 
the luminosity function which would be measured in a galaxy survey with a 
given limiting surface brightness. We show that current measurements of the 
galaxy luminosity function may be missing more than half of all galaxies, 
and an even larger fraction of faint galaxies, given the correlation between mass 
and surface brightness. The likely underestimate of the luminosity density is 
also expected to be large. We discuss how this affects observations of the "faint 
blue galaxy" population. 

We also investigate the dynamics of galaxies as a function of surface 
brightness. We show that, in the absence of any systematic change in the ratio 
of disk mass to disk luminosity, galaxies of all surface brightnesses should lie on 
the same Tully-Fisher (1977) relation. Our models also show systematic changes 
in the shape of the rotation curve as a function of angular momentum, which 
leads to low surface brightness galaxies having slowly rising rotation curves. 
Furthermore, because high angular momentum LSB disks have their baryonic 
mass spread over a larger area than normal galaxies of comparable mass, LSB 
disks contribute very little to the observed dynamics of the galaxy. Thus, LSBs 
provide a very effective tracer of the shape and mass profile of the dark matter 
halo, out to proportionally larger radii than is possible to observe with normal 
galaxy rotation curves. 



1. Introduction 

Disk galaxies, including low surface brightness galaxies (LSBs), show remarkable 
regularities: most disk galaxies have asymptotically flat rotation curves (see Bahcall and 
Casertano 1985, Burstein & Rubin 1985, Persic et al. 1995, de Blok et al. 1996, and 
the references therein); the Tully-Fisher relation between circular velocity and luminosity 
appears to hold over a wide range of both mass and surface brightness (Strauss and Willick 
1995, Zwaan et al. 1995, Sprayberry 1995); and the light profile of galactic disks are well fit 
by an exponential over several disk scale lengths (de Jong 1995) for high surface brightness 
spirals, as well as for non-dwarf LSBs (Davies, Phillips, & Disney 1990, McGaugh & 
Bothun 1994, James 1994, McGaugh et al 1995, Dalcanton et al. 1997, DeJong 1995). The 
consistency of these properties suggests that the basic structure of all disk galaxies is a 
general property of galaxy formation, either tied to regulated feedback during formation, or, 
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as we shall consider in this paper, to primordial properties of protogalaxies, namely their 
mass and angular momentum. The fact that galaxies' rotation curves are smooth in both 
the disk-dominated inner regions and the halo-dominated outer regions strongly suggests 
that the structure of the inner halo is coupled to the formation of the disk. 

Past disk formation models have been able to reproduce many of the universal 
properties of disk galaxies (Fall & Efstathiou 1980, Gunn 1982, van der Kruit 1987). 
The models were motivated by Mestel's (1963) observation that the angular momentum 
distribution of the galactic disk is very similar to that of a sphere in solid body rotation. 
Thus, the models assumed that the collapse of a uniformly rotating gaseous protogalaxy was 
a good model for the formation of the disk. In general, the models produced realistic disk 
surface density profiles, provided that the models included the further assumptions that 1) 
there was very little angular momentum transport while the gas is collapsing, 2) the angular 
momentum was responsible for eventually halting the collapse, and 3) the density profile 
of the dark matter halo was such that the resulting galaxy will have a flat rotation curve. 
However, because of the third assumption, the resulting mass distribution of the disk+halo 
system was not necessarily self consistent; because the halo did not respond to the mass of 
the collapsed disk, the final disk was not necessarily in gravitational equilibrium. 

One limitation of past models of disk formation is that they have only attempted 
to reproduce the properties of "normal", high surface brightness disks such as our own. 
However, a growing body of literature on low surface brightness galaxies (LSBs) has 
revealed that there is a substantial population of disk galaxies whose photometric properties 
fall well outside what is usually classified as "normal" . Any viable formation scenario must 
explain the tremendous range of disk galaxy surface brightnesses and scale lengths, while 
simultaneously explaining the tremendous similarities among disk galaxies throughout this 
range. 

In this paper we develop a gravitationally self-consistent scenario for the formation 
of disks It is remarkably effective at explaining both the structural and dynamical 
regularities among disks, as well as for motivating the enormous range of galaxy properties 
which are observed. We use the model to calculate directly observable properties of 
individual galaxies (light profiles [§ |2.1.1|| , rotation curves [§ [2.1.2|| , etc.). We also use it to 



compute the expected distribution of galaxies as a function of surface brightness and disk 
scale length; this allows us to easily calculate the degree to which selection effects shape the 
observed properties of galaxies [§E|. 
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2. The Formation of Disks 

The outline of the scenario which we will use for disk formation is as follows: 

(1) In the early universe, tidal torquing spins up both the dark matter and the baryonic 
matter in some region of space, which we will approximate as an isolated sphere of radius a 
in solid body rotation (i.e. a rotating, spherical top-hat overdensity) . While the region has 
enough of a mass overdensity that it will eventually collapse into a galaxy, the overdensity 
is small at early times (5 cx t 2 ^ 3 ) and the mass density can be treated as uniform to first 
order. We also assume that the gas and the dark matter are uniformly mixed. Because 
inhomogeneities in the density field are sufficient to provide tidal torques (Peebles 1969), 
we do not require that large numbers of galaxies must have condensed before tidal torquing 
can be effective. 

(2) As time passes, both the dark matter and baryonic matter begin to collapse in the 
overdense region. Initially, when the densities are low, the baryons cannot cool, and the 
pressureless collapse proceeds identically for both components. 

(3) Because the dark matter cannot dissipate energy, its collapse halts when the system 
virializes. Numerical simulations suggest that the resulting halo has a generic density profile 
(Warren et al. 1992; Dubinski and Carlberg 1993; Navarro et al. 1996): inside the core 
radius, the density profile rises roughly as r" 1 , while outside the core radius, the density 
profile falls as r -7 , with 7 between 3 and 4. We will approximate the halo with a Hernquist 
density profile, whose asymptotic radial profile scales as r~ 4 (see eqn [I] below). The size of 
the core radius depends only on the mass of the halo and the mean density of the universe. 
For a flat universe (which we will assume for simplicity for the duration of this paper), the 
growth of the halo is self-similar, and the only change in the halo with time is the increase 
in the core radius due to its increasing mass. 

(4) At some point during the joint collapse of the dark matter and baryons, the density 
over some region becomes high enough that the baryons begin to cool and decouple from the 
dark matter. As they cool, their collapse accelerates, and the baryons begin to concentrate 
within the dark matter. This condensation progressively increases the baryonic fraction 
within the inner parts of the halo. Because the baryons have non-zero angular momentum 
they cannot collapse all the way to the center, and instead settle into a rapidly rotating 
disk. The final mass distribution in the disk is determined both by the initial distribution 
of specific angular momentum and by the final rotation curve of the collapsed disk+halo 
system. 

(5) After the onset of cooling, the collapsing baryons further condense the dark matter 
halo by increasing the mass density in the inner parts of the halo, where the baryonic 
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mass fraction has increased the most. This modifies the Hernquist potential which the halo 
had at the time when cooling began. The condensation of the halo is assumed to proceed 
roughly adiabatically. 

The calculation outlined above, which we perform explicitly in the following section, 
expands upon the formalisms of Fall & Efstathiou (1980), Gunn (1982), and van der Kruit 
(1987) by dropping the assumption of a static halo and allowing the disk halo to be pulled 
in by the collapsing baryons, while retaining the initial condition of a uniformly rotating 
protogalaxy and the assumption that the collapse time is much shorter than the angular 
momentum transport time. 

An essential ingredient in our calculations is the assumption that angular momentum 
transport is negligible during the formation of disk galaxies. We have several reasons 
to believe that angular momentum transport is generally small, whether the formation 
is smooth, as we have assumed, or hierarchical, taking place through merging of low 
mass subunits. As we discuss at length in §^|, the observed angular momenta of disks is 
comparable to what is expected for galaxy halos, suggesting that there has not been a 
dramatic drop in the angular momentum of the gas during formation. Furthermore, in 
numerical simulations within which the gas fails to conserve angular momentum during 
galaxy formation (e.g. Navarro Sz Steinmetz 1996), the disks which form are too compact 
and have too little angular momentum. In other words, angular momentum transfer 
between the gas and the dark matter halo leads to the formation of galaxies whose density 
profiles differ strongly from those of real galaxies. This empirical argument suggests that 
angular momentum transport between the gas and the halo is generally not important for 
the formation of disk galaxies. 

Angular momentum transport during collapse is most likely to be negligible for the 
formation of low surface density disks, which as we shall show in §[3], are the dominant 
form of disk galaxy; angular momentum transport is likely to be most effective in 
baryon-dominated systems such as ellipticals and bulges and least effective in low surface 
brightness systems. A pressure supported dark halo tends to suppress the formation of 
both bars (Ostriker and Peebles 1973) and spiral arms (Athanassoula 1984), two effective 
means of transporting angular momentum. In low surface density, high angular momentum 
systems, the dark matter dominates the potential, so bars and spiral arms are unlikely to 
form. On the other hand, in low angular momentum systems, dissipation concentrates the 
baryons to high densities, where they are subject to non-axisymmetric instabilities that 
will transport angular momentum. Thus, our basic assumptions are most likely to be valid 
for LSBs and are almost certainly not correct for ellipticals and the bulges of high surface 
brightness galaxies. Finally, we expect very little angular momentum transport within the 
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disk after its formation. The angular momentum transport time due to spiral density waves 
is quite long, ~ 5 — 10 Gyr (Zhang 1996, Gnedin et al. 1995), and thus in the absence 
of bar formation or other large scale instabilities, disks cannot redistribute their angular 
momentum efficiently once formed. 

The other major assumption about angular momentum is that the baryons have an 
initial angular momentum distribution similar to that of a uniformly rotating sphere. As 
in past disk formation models, our assumption is motivated by the similarity between the 
observed disk angular momentum distribution and the angular momentum distribution of a 
uniformly rotating sphere (Metsel 1963). Furthermore, in a simple picture of tidal torquing, 
most of the torqueing should occur near maximum expansion and will be due to an external 
quadrupole. This external torque will produce uniform rotation. Our assumption is also 
testable by comparision to numerical simulations, although current simulations have not 
addressed this problem explicitly. 

Finally, in our model we appear to have assumed that all of the baryons collapse 
simultaneously, which is physically unreasonable. The baryons will begin to collapse beyond 
the dark matter when they begin to dissipate energy. The onset of cooling, however, is a 
density dependent phenomena, and the baryons in the densest, inner region of the halo will 
begin to collapse well before the outer region. However, because we are assuming that the 
system is spherically symmetric to first order, at any given time both the baryons and dark 
matter outside of the cooling region will be largely unaffected by how the matter inside 
the cooling region rearranges itself. Thus, while the arguments below are dependent on 
the initial and final distribution of the baryons, they are independent of the details of how 
the collapse proceeds to its final state, as long as the process is spherically symmetric and 
adiabatic. 



2.1. Simple Model 

In this section we will use the above scenario to calculate the mass profile of the disk 
and halo after collapse. The evolution of the mass profile can be traced through three 
distinct stages, each identified with a distinct length scale: 1) the uniform sphere before 
collapse, with radius a; 2) the Hernquist profile with core radius r , at the moment when 
the majority of the baryons begin to decouple from the dark matter; and 3) the final 
baryonic disk, characterized by an exponential disk scale length a. First we will use energy 
conservation to relate the core radius of the Hernquist profile halo to the radius of the 
initial protogalaxy. Next we will use adiabatic dragging of the halo and the initial angular 
momentum distribution of the disk to solve for the final mass profile and rotation curve 
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of the resulting galaxy. Finally, we will find the scale length, central surface brightness, 
and circular velocity of the resulting baryonic disk as a function of baryonic mass fraction, 
initial angular momentum, and the total galaxy mass. 

The core radius r of the Hernquist profile can be easily related to the initial 
protogalaxy radius a, using energy conservation. Letting M tot be the total mass within the 
initial sphere of radius a, and F be the fraction of mass within a which is baryonic, the 
potential energy of the dark matter in the initial sphere is —3(1 — F)GM^ ot /5a. This must 
be equal to the potential energy of the dark matter when the galaxy has partially collapsed 
into a Hernquist mass profile, defined as: 



where Mhaio (r) is the halo mass within radius r. Its potential energy within a radius r is 



Fixing the mass within r = a to be constant, = M tot [(a + r )/a] 2 . Assuming that 
the energy of the initial uniform sphere is entirely gravitational, that the energy of the halo 
within a is conserved during collapse, and that the energy of the collapsed Hernquist halo is 
roughly consistent with virialization (E rs |$| /2), one finds that a = 3.2r . If virialization 
is not complete, the ratio of a/r$ will be somewhat smaller, but will not significantly change 
any of our results. The relationship between a and tq, as well as all of the other relations 
derived in the section, are rederived in the Appendix for a shallower halo potential suggested 
by Navarro et al. (1996). With the shallower potential, a = 7.1r . For a flat universe, the 
evolution of a dark matter halo is a scale free process, thus, the density at the core radius 
should be a constant multiple of the mean density of the universe: r = C~ l M^p 1 ^ 3 , 
where C ~ 70 (Navarro et al. 1996) and p is the mean density of the universe. This allows 
the mass scale of the galaxy to be related to the initial size of the perturbation. 

The next step is to calculate the mass distribution of the final disk, by considering 
the angular momentum distribution of the baryons and the adiabatic dragging of the dark 
matter by the collapsing disk. We have assumed that the baryons are evenly mixed with 
the dark matter before the onset of cooling, and that they have the distribution of specific 
angular momentum characteristic of disks today, which is well represented by a sphere in 
solid body rotation (Mestel 1963): 




(1) 
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M gas (< j) = FM tot [l - (1 - j/j max f 2 ] , (3) 

where M gas {< j) is the baryonic mass of the galaxy which has specific angular momentum 
less than j, and where j max = Qa> 2 for a sphere in solid body rotation with angular rotation 
velocity fl The maximum angular momentum can also be characterized by the spin angular 
momentum parameter A = Jtot\E\ l ^ 2 G~ l M~ 5 / 2 (Peebles 1969); for our initial conditions (i.e. 
an unvirialized sphere in solid body rotation) J tot = 2M tot j max /5 and E = —3GM 2 ot /5a, 

1 /2 

and therefore j max = § (§) \\fGM tot 



a. 



We can now solve for the final mass distribution of the disk+halo system. To 
characterize the radial distribution of halo mass, we define m^(r) = M/ l(l / (r)/(l — F)M to t 
to be the fraction of the halo mass enclosed within a radius r. We will use rrih as a 
radial parameter for the rest of this calculation, as well as a similar parameter related 
to the enclosed disk mass, m d {r) = M disk {r) / F M tot . (We will henceforth drop the 
explicit dependence of rrih and m d upon r for notational convenience.) Note that the 
radial dependence of these parameters, r(m^) (or r(m d )), will change during collapse, 
and thus while always defines a surface bounding a certain fraction of the halo 
mass, the radius of that surface is not constant. After collapse, we will let M(m^) 
(= M to t [(1 — F)mh + Fmd(mh))) be the total mass enclosed within r{mh)- 

Assuming that conservation of angular momentum is responsible for finally halting the 
collapse of the disk, the specific angular momentum of gas which halts at radius r(md) is 
j( m d) — \jGM(md)r(m ( i). All gas which has specific angular momentum less than j(m^) 
will wind up interior to r(m^), and thus, using eq. |3|, 



f Mtot 

3/2 



GM(m d )r(m d ) , 

(5) 



Now, to calculate the final mass profile of the halo, the adiabatic invariance of the 
angular action Ig (= / vg ■ r dQ) must be used (see Binney & Tremaine 1987, Chapter 
3). Because the baryons collapse by a large factor, the time scale for their collapse is 
much larger than the orbital time of the dark matter at a given radius; therefore, the 
collapse is roughly adiabatic, and the orbits of dark matter particles will adjust to the 
changing potential such that If = GM(r)r is constant (Blumenthal et al. 1985, Flores et 
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al. 1994). Using the inverted form of the Hernquist potential for the initial conditions, 
r i( m h) = r o m // 2 /(l + r o/ a ~ m // 2 )> after collapse the angular action is 



GM(m h )r(m h ) = Gm h M tot ri{m h ) (6) 

3/2 

777 

= GM tot r — —h- - 2 . (7) 

1 + {ro/a) - m h ; 

Because M(mh)r(mh) = M(rnd(rnh))r(md(rnh)), eqn |G] may be substituted into eqn | to 
solve for m^rn^), which then yields 

/ t 3/4 \ 3 / 2 

m d (m h ) = 1 - 1 - , , (8) 

V ^l + (r /a)-ml /2 J 

where £ = V GM tot r / j max , or £ = 0.173/A. Note that for small enough values of A, there is 
some maximum for which eqn. |8| is well defined; this cutoff corresponds to the maximum 
extent of the collapsed disk. For small values of lambda, the disk collapses to well inside 
the final extent of the halo, and beyond this extent, should be fixed at 1. 

The angular action may also be used to derive the final radius enclosing rrih of the halo 
mass: 



1 h > GM(m h ) 1 l + Wa )_„f \(l-F)m h + Fm i (m h ) ] 



Coupled with 



M(m h ) = M tot [(1 - F)m h + Fm d (m h )} , (10) 

equations §, ^|, and [lOj are parametric equations which fully specify the mass profile of the 
the collapsed galaxy. The final mass profile of the disk and halo therefore depends only 
on the baryonic mass fraction F, the density of the universe po> an d the mass M tot and 
spin angular momentum A of the protogalaxy. Note that the only place where our choice 
of initial conditions (i.e. the uniform sphere) plays a significant role is in determining the 
numerical factor relating £ to A. 



-10- 



The above equations can be solved to give both the surface density and the rotation 
curve for the final disk: 

yj( r ( m \\ _ FM tot dm d dm h 
2nr(mh) drrih dr 

and v1{r-{m h )) = GM{m^)/r{mh). 

As we will show in the following sections, for a 1O 12 M halo with F = 0.05 and 
A = 0.06, the above model implies r = 35kpc, L = 1.7 x 1O 1O L , t> = 250 km/s, a = 3.4 
kpc, and /j,q = 21.2 mag/arcsec 2 , all reasonable parameters for a massive disk galaxy 
(assuming a mass-to- light ratio of 3 M / L for the baryons in the disk). 



2.1.1. Surface Density 

The calculations above solve for the mass of both the dark matter and the baryons as 
a function of radius, and thus they can be used to calculate the surface density profile of 
the baryonic disk which results from the collapse. Figures |l|a and |I]b show the disk surface 
density profile which results from eqn. [11] for different values of the dimensionless spin 
angular momentum A and the galaxy mass M tot . All of the disk profiles produced by this 
gravitationally self-consistent model are roughly exponential over many scale lengths, with 
a scale length a that depends upon A, M to t, F, and tq. After fitting the disk profile between 
0.5 and 2.5 scale lengths for a wide range of parameters, we derive this approximate fitting 
formula for the scale length: 



16F 2 + 8F + 1 + (0.002/F) . 

(12) 

which is good to ±5% over 0.03 < A < 0.18, 0.02 < F < 0.15, and 1 kpc < r < 120 kpc; 
for larger values of A, the formula overestimates the scale length by ~ 10% at A = 0.2 and 
~ 30% at A = 0.4. There is an apparent cut-off in the disks at large radii in Figure |l] which 
is due to using a sharp edged sphere in the initial conditions 

Using the fact that the total disk mass is the baryonic fraction F times the total galaxy 
mass Mtot, the approximate central surface brightness is 



a(M to t,X,F) ~2.55r (M, 



tot J 



jL-fa.r 



2 + A 



-i 



(0.084/A) 5 
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So(M - A ' F) = ^ZIfyt (13) 



13ttT 



(1 + 8F + 16F 2 + ^) 2 



M 1/3 

m tot 



2 + A- x + 



0.083 
A 



2 -| 2+6F 



(14) 



where T = [Mdiskl Afofc)j note that T is not a dynamical mass-to-light ratio, but a measure 
of the efficiency of turning gas into starlight over the lifetime of the disk. 

To first order, S oc F X~^ 2+6F \ suggesting that low mass and/or high angular 
momentum protogalaxies naturally form low surface brightness galaxies. Physically, because 
of angular momentum conservation, a more rapidly rotating protogalaxy will be unable 
to collapse as far inwards as a less rapidly rotating one. Thus, the galaxy with a larger 
specific angular momentum will spread its baryonic mass over a larger disk, leading to 
lower mass surface densities (as seen in eqns. [12| and [14] ) . There is observational support 
for this picture from Zwaan et al. (1995) who find that LSBs have systematically larger 
scale lengths at a given circular velocity than their high surface brightness counterparts; 
given that the specific angular momentum must scale like aV c (oo) to first order, LSBs must 
have higher angular momenta than normal galaxies of the same mass. Most importantly, 



Equation |14] also suggests that any spread in either the masses or the angular momenta of 
protogalaxies translates into a spread in the resulting distribution of disk surface densities. 
We will explore this effect in §||. 

One feature of the models is the presence of a central cusp in the surface density 
distribution (Figure []]). The mass in the cusp is only a small fraction of the total disk mass 
(~ 6%), and is largely independent of the choice of the initial mass profile of the halo (see 
Appendix A). The cusp is a consequence of our assumption that the angular momentum 
distribution is well-modelled by a uniformly rotating sphere (equation |3]). The center of the 
disk contains the very low angular momentum material from the initial protogalaxy, and 
for a uniformly rotating sphere, M(< j) oc j for j << j max , implying S oc r _1 in the central 
regions where the halo contributes very little to the total mass. In contrast, the surface 
brightness distribution of a pure exponential disk is constant for very small radii, and thus 
the surface brightness profiles that result from our models must necessarily be steeper than 
exponential in the center. Because the cusp is an artifact of our simple initial conditions, 
we do not ascribe any profound importance to the bulge-like centers of the model's light 
profiles. Any variation in the shape of the initial perturbation, the alignment between the 
angular momentum and the principle axes of the perturbation, or the degree of angular 
momentum conservation will lead to variations in the size and shape of the cusp, and thus 
suggest natural mechanisms for producing variations in the bulge-to-disk ratio. 
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Fig. I.- 

Resulting surface density profiles [a,b], rotation curves [c,d], and disk mass profiles [e,f] for 
disk galaxy models with various angular momentum (left column; A = 0.03 — 0.18, linearly 
spaced, M to t = 10 12 M ) and mass (right column; M to t = 10 10 — 10 12 M , logarithmically spaced, 
A = 0.06). In all panels, the heavy line is a fiducial model with A = 0.06 and M to t = 10 12 M ; 
we have assumed F = 0.05 and H$ = 50km/s/Mpc for all the models. The top row [a-b] shows 
the surface density as a function of radius, with the right hand axis giving the apparent surface 
brightness /iq assuming a disk mass-to- light ratio of T = 3Mq/Lq. The horizontal dotted 
lines bracket the Freeman (1970) surface brightness. The solid lines are the models, the dotted 
diagonal lines are exponential fits to the surface brightness distribution, and the solid dots are 
at one exponential scale length (r = a). The models produce roughly exponential profiles over 
several scale lengths. Panel [a] shows how increasing the angular momentum of the protogalaxy 
decreases the extrapolated central surface brightness and increases the exponential scale length 
of the resulting disk, while panel [b] shows how increasing the mass of the protogalaxy increases 
both the central surface brightness and the scale length of the resulting disk. The middle row 
[c-d] shows the rotation curves which result from the models. Panel [c] shows that changing the 
angular momentum changes the shape of the rotation curve, such that high angular momentum 
models rise more slowly, but keeps the same asymptotic circular velocity. Panel [d] shows that 
decreasing the mass decreases the asymptotic circular velocity; the shape of the rotation curve 
does not change however, and simply rescales with the exponential scale length and mass. The 
bottom row [e-f] shows the baryonic-to-total mass fraction within a given radius (solid lines), as 
well as the fraction of the baryonic mass contained within a given radius (dotted lines). These 
plots show the degree to which the dynamics of the galaxy are strongly affected by the baryons. 
Panel [e] shows that in low angular momentum models (A < 0.04), almost all of the baryonic 
mass lies within a radius where baryons make up more than 40% of the total mass, whereas in 
high angular momentum models (A > 0.09), almost none of the baryonic mass falls in a region 
where baryons make up more than 40% of the total mass. The models also suggest that in the 
high angular momentum, low central surface brightness models, the halo dominates the dynamics 
over almost the entire disk. Panel [f] shows that for changing mass, the baryonic mass fraction is 
self-similar, scaling with the exponential scale length. 
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2.1.2. Rotation Curves 

Equations || and [H] can readily be used to calculate the circular velocity as a function 
of radius. Figures |I|c and |T]d show the resulting rotation curves for the same disk plus halo 
models shown in Figures |l|a and |l]b respectively, for different values of the spin angular 
momentum A and the total mass M to t- For simplicity, we have included only the monopole 
term from the disk and not thus included the effects of disk flattening. The rotation curves 
confirm Blumenthal et al.'s (1985) and Flores et al.'s (1995) conclusion that adiabatic 
contraction of the halo leads naturally to asymptotically flat rotation curves, and, as 
previous authors have found, that for a wide range of parameters, the rotation curves of the 
models look much like those of real galaxies. 

Figure |l][c&d] demonstrates how changes in mass and angular momentum affect the 
shape of the rotation curve V c (r). First, the shape of the rotation curve remains constant 
with changes in mass, when scaled by the exponential scale length, for fixed spin angular 
momentum A. At fixed angular momentum, the baryons collapse by the same factor within 
the self-similar dark matter halos, leading to identical fractions of dark and luminous matter 
within a given dimensionless radius r/ro (or r/a), and thus to self-similar rotation curves. 
However, if there is any deviation from perfect self-similarity in the shape of dark matter 
halos, then there will be accomanying deviations in the self-similarity of the rotation curves. 
Violations of our assumption of perfect self-similarity will lead to either increased scatter in 
the shape of rotation curves at fixed angular momentum, or to systematic variations in the 
shape of rotation curve with changing mass. 

In contrast, changes in the spin angular momentum A lead directly to variations in the 
shape of the rotation curve. Increasing the spin angular momentum of a galaxy decreases 
the collapse factor of the baryons (i.e. rj/ry), leaving them at a larger dimensionless radius. 
Thus, an increase in the spin angular momentum reduces the fraction of baryonic mass 
within any given dimensionless radius, leading to a rotation curve whose shape is dominated 
more by the dark matter distribution than the baryonic matter. There are two strong 
observational signals of these high angular momentum rotation curves. First, high angular 
momentum galaxies will have a larger dynamical mass-to-light ratio, due to the increased 
dominance of dark matter over baryons at every radius (as can be seen in the mass profiles 
of Figure |l|[e&f]). Second, the rotation curves will rise gradually with radius, instead of 
steeply as is usually observed for more concentrated low angular momentum disks. 

As shown in § |2.1.1| , the models suggest that such high angular momentum systems 
should to be found preferentially among the population of low surface brightness galaxies 
(eqn [14]). Indeed, the rotation curves of LSBs do reveal the same signatures as the 
high-angular momentum models. Edge-on (Goad & Roberts 1981) and face-on LSBs (de 
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Blok et al 1996) show a systematically slower rise in their rotation curves than do their 
high surface brightness counterparts. De Blok et al (1996) also find that the dynamical 
mass-to-light ratio increases systematically with surface brightness, in concordance with the 
model. 

The rotation curves in Figure 0[c&d] and the mass profiles in Figure |I|[e&f], suggest 
that low surface brightness disks are extremely good tracers of the properties of dark 
matter halos. First, because low surface density disks contribute less to the dynamics of the 
galaxy within any given radius, the resulting rotation curve is much closer to what would 
be expected for a massless disk. Second, LSB disks will tend to reflect any deviations from 
axisymmetry within the halo. Because the disks do not have enough mass to round out the 
potential at a given radius, the disks should be more asymmetric in general; however, if 
LSBs are high angular momentum systems, they may preferentially form in halos with large 
quadrapole moments, which will lead to a further increase in the apparent disk asymmetry. 
Finally, LSBs have larger scale lengths than high surface brightness galaxies of the same 
mass. This allows the halo to be traced to a much larger radius, allowing the mass profile 
of the halo to be measured at extremely large distances. 

We may make a detailed comparison between the shapes of observed rotation curves 
and the rotation curves produced by the models. Figure || shows the logi?-logM(_R) 
relations for the Burstein & Rubin (1985) "mass types" (derived from observed rotation 
curves) and the identical relation for models spanning a moderate range of spin angular 
momenta. (All of the curves have been rescaled by an arbitrary factor R m and M m to 
match at their point of maximum inflection (as in Burstein & Rubin 1985), and the model 
curves have been shifted upwards from the observations to make the comparison easier.) 
Figure shows that, one, both the models and the observations are remarkably similar, 
with the rotation curves being well approximated a broken power-law, and two, that it is 
the faster rising rotation curves which have less curvature. In the models, these changes in 
shape are produced by reasonable variations in the spin angular momentum. Furthermore, 
because changes in mass do not produce changes in the shape of the rotation curve, there 
should not be any strong correlation between galaxy mass and the distribution of Burstein 
& Rubin "mass types" , unless the dark matter halos are not self-similar with mass, as we 
have assumed. There is a tendency for the model mass distributions to be slightly shallower 
at small radii than the observed distributions, which probably reflects either a deviation 
from the assumed initial Hernquist halo mass profile or in the initial distribution of angular 
momentum. Given the simplicity of the initial assumptions, this is less notable than the 
general degree of agreement. 
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Fig. 2. — Comparison between the shapes of the rotation curves resulting from models 
with A = 0.03 (light solid), A = 0.05 (medium solid), and A = 0.03 (heavy solid), and the 
Burstein & Rubin (1985) "mass types" I (light dashed), II (medium dashed), and III (heavy 
dashed). The curves have been rescaled in radius and mass to match at their points of 
maximum curvature, as in Burstein & Rubin (1985). Changes in angular momentum are 
clearly capable of producing variations in shape comparable to those found observationally 
by Burstein & Rubin. The models rise less steeply in the center and more slowly in the 
outer regions than do the Burstein & Rubin mass types. This probably reflects a difference 
between reality and our assumed halo profile or angular momentum distribution, as well as 
possibly our neglect of bulges. 



-17- 



Recent work by Persic, Salucci, & Stel (1996; see also Persic & Salucci 1996) has 
parameterized the shape of galaxy rotation curves as a function of mass to produce a 
"Universal Rotation Curve" (URC). However, as our modelling shows, such attempts 
may perhaps reveal more about the selection criteria of the sample of galaxies used to 
generate the URC than about the rotation curves themselves. First, because mass and 
angular momentum are fundamentally linked with surface brightness and exponential 
scale length, any selection bias which affects the distribution of surface brightnesses and 
scale lengths present in a galaxy sample will critically affect any URC derived from the 
sample; any selection effect which changes the distribution of surface brightnesses and 
scale lengths in a sample can change the apparent distribution of disk angular momenta, 
and thus the distribution of rotation curve shapes. For example, intrinsically large scale 
lengths will be overrepresented in angular diameter limited surveys of field galaxies, while 
an angular diameter limited surveys of cluster galaxies will have a larger fraction of 
galaxies with small scale lengths. Thus the two surveys will have different distributions of 
galaxy masses and angular momenta, giving different apparent relations between mass and 
rotation curve shape. Furthermore, there are unavoidable biases against finding low surface 
brightness galaxies, which will tend to reduce the contribution of low mass or high angular 
momentum galaxies to the derived URC. To avoid drawing any erroneous conclusions, any 
parameterization of rotation curve shapes must include a detailed discussion of the selection 
criteria that went into choosing the galaxies whose rotation curves were measured. 



In addition to understanding the shape of galaxy rotation curves, our models may 
be used to understand the behavior of the Tully-Fisher relation as a function of surface 
brightness. Towards that end, fitting the circular velocity at three disk scale lengths, we 



2.1.3. Tully-Fisher Relation 



find 



V c (r = 3a; r , A, F) — 8 km/ s 



r (M tot ) 
kpc 



[a(X)F + b(X)} 



(15) 



where 



a(A) = 1 + 35 



(0.015/A) 2 



(16) 



1 + (0.015/A) 2 



and 
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6(A) = 0.22 + 0.68 



(A/0.015) 
1 + (A/0.015) 



(17) 



which is accurate to ±5% for 0.01 < A < 0.14, 0.01 < F < 0.15, and lkpc < r < 120 kpc. 



The lambda dependence (i.e. the term in brackets) in equation [15] contributes only a 
±20% variation in V c for A > 0.025. For very small angular momentum models, however, 
the concentration of the disk is substantial enough that the mass of the compact disk 
contributes significantly to the rotation curve, giving a factor of two variation in V c for 
A < 0.025. However, these systems are the least likely to survive as pure disk galaxies; 
these low angular momentum disks are disk dominated in their centers and are likely to be 
unstable to bar formation. The bar would transfer angular momentum outwards to the halo 
(Weinberg and Hernquist 1992), which redistributes angular momentum within the galaxy 
and alters its structure. The low angular momentum systems are also plausible candidates 
for elliptical galaxies. 



lU, =19-25 T=3 M S /L Q F=0.05 H = 50 km/s/Mpc 




V (km/s) 



Fig. 3.- 

Comparison between the Tully-Fisher relation predicted by equation [l5]for /j,q = 19 — 25 mag/arcsec 2 (solid 
lines, lighter corresponding to fainter surface brightness) and the observed Tully-Fisher relation (Her) for 
clusters (dashed line; Willick et al. 95) and the field (dotted line; Willick et al. 96). Note that, for the 
assumption of a constant mass-to-light ratio for galactic disks, there should be no shift in the mean Tully- 
Fisher relation as a function of surface brightness. We have assumed T = 3M Q /L Q , F = 0.05, and 
Hq = 50km/s/Mpc for the models. 
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Equation |TJ above can be used to understand the Tully-Fisher relationship (Tully 
& Fisher 1977; see Strauss & Willick 1995 for a review) for disk galaxies, for all surface 
brightnesses. Because the scale length of the disk is proportional to r (eqn. |T2D , and the 
circular velocity at fixed scale length is proportional to ^/M(r)/r , giving L cx V^, with 
a numerical constant that depends upon F, Y, and A. This slope is consistent with the 
observed infrared Tully-Fisher relation (Strauss and Willick 1995). Furthermore, because of 
the weak dependence of V c upon A, equation [15| predicts very little offset in the Tully-Fisher 
relation with changes in A, and thus predicts that the Tully-Fisher relations for both high 
and low surface brightness galaxies should be nearly coincident. We show this in Figure |3], 
where we plot the relationship between circular velocity and absolute magnitude predicted 
by equation [15] for a wide range of disk surface brightnesses. Superimposed are the ±ler 
limits for the Tully-Fisher relation derived by Willick et al (1995, 1996) for normal cluster 
and field galaxies. The predicted Tully-Fisher relations fall well within the observed scatter 
over a factor of 100 in surface brightness and a factor of 1000 in mass. 

The independence of the predicted Tully-Fisher relation and disk surface density is not 
surprising given that the relation primarily results from the properties of the dark matter 
halo, if one reasonably assumes that the luminosity of a galaxy is proportional to the total 
mass, and that the measured line width reflects the asymptotic circular velocity at large 
radii, where the mass is dominated by the halo; the relationship between luminosity and 
line width should not depend strongly on how the baryons are distributed within the halo. 
However, the predicted relations do assume a single disk mass-to-light ratio T for all masses 
and all disk surface densities; any variation in T should translate into increased scatter in 
the relation, and any systematic variation in T with the total mass or surface density of 
the disk should translate into changes in the Tully-Fisher slope. Therefore, the remarkable 
agreement between the predicted Tully-Fisher relation in Figure |] and the measured one 
implies that there is very little evidence for a systematic variation of the disk mass-to-light 
ratio with mass. On the other hand, while Sprayberry et al. (1995) and Zwaan et al. 
(1995) both find that LSB galaxies fall on the same Tully-Fisher relation defined by normal 
galaxies, they also find increased scatter in the relationship, roughly by factors of two. This 
suggests that there is more variation in the disk mass-to-light ratio for low surface density 
disks than for high surface density disks; given the more tenuous nature of these disks, we 
do not find this conclusion to be unreasonable. However, given that the mean Tully-Fisher 
relation is indistinguishable for both high and low surface brightness galaxies, the mean 
disk mass-to-light ratio must be similar across the range of observed surface brightnesses. 



2.2. Pressure Support 



In the preceding calculation of the galaxy formation, we assumed that the baryons 
will collapse along with the dark matter halos. However the baryonic gas experiences 
pressure forces in addition to gravitational forces. For small masses, internal pressure may 
support the gas against collapse, in spite of the inward gravitational pull of the collapsed 
dark matter halo (Jeans 1929). Because of the link between low masses and low surface 
brightnesses (eqn [14]), there is some limiting surface density for which we expect our 
assumption of complete baryonic collapse to become invalid. 

To contrain the regime of galaxy properties for which pressure support can affect 
the formation of a disk, we consider the simpler question: what is the range of halo 
properties for which the gas is bound to the dark matter halo? The gas will be bound if the 
gravitational binding energy is greater than its thermal energy. Assuming that the gas is 
fully ionized, it will be gravitationally bound to the dark halo if 



^2 > kT, (1* 
2 ' v 

where V c is the circular velocity of the dark matter halo, m p is the proton mass, and T 
is the temperature of the gas. Applying this relation at the halo half mass radius, where 
v c = JGM/vq = V GCM 1 / 3 p ^ 6 implies a minimum halo mass for collapsed objects: 



M min = po ~ ~ = 4 x 1O 7 M T 4 3/2 (19) 

\Gp m p C J 

where T4 = T/10 4 K. Thus, the Jeans criterion limits the formation of low mass galaxies, 
rather than low surface brightness galaxies, and we can treat the above calculation as valid 
for all protogalaxies with M tot > M min . 



3. The Number Density of Disk Galaxies 

The equations of §2.1.1| provide a means to relate the primordial quantities of mass and 
angular momentum to the observable properties of disk galaxies, namely central surface 
brightness and scale length. Theoretical models and numerical simulations provide a way 
to predict the distribution of mass and angular momentum for galaxy halos, and thus we 
have the means to predict the observed distribution of disk central surface brightnesses 
and scale lengths (a). Furthermore, because most selection effects in galaxy catalogs can 
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be modelled as detection efficiencies as a function of surface brightness and/or scale length, 
the resulting joint distribution of /io and a can be convolved with a model of a survey's 
detection efficiency to predict the observed distribution of either /i and a, or equivalently, 
surface brightness and luminosity. 

To begin, we must first assume an intrinsic distribution for spin angular momenta, 
p(A). This problem has been approached both analytically and numerically by a number 
of groups (Barnes & Efstathiou 1987, Ryden 1988, Warren et al. 1992, Eisenstein & Loeb 
1995, Catelan & Theuns 1996a, Catelan & Theuns 1996b), all of which find the same 
general properties. To first order, the distribution of spin angular momentum of collapsed 
dark matter halos is well approximated by a log normal distribution: 



Most work finds that this distribution peaks around (A) pa 0.05, with a width in the 
log of a\ ~ 1. For the purposes of this paper, we adopt the specific values of a\ ~ 0.7, and 
(A) ~ 0.06, which provide a good fit to the results of Warren et al. (1992). 

For simplicity, we will assume that p(X) is identical for all sizes of the initial density 
perturbation (i.e. angular momentum is independent of mass). Faber (1982) has argued 
that the spin angular momentum of a protogalaxy should be independent of the amplitude 
of the initial overdensity, based upon scaling laws. However, there is some evidence from 
analytic calculations of p(X) in the peak-height formalism which suggests that the mean 
spin angular momentum, (A), is a function of the height of the peak in the initial density 
field which forms a given protogalaxy (Eisenstein & Loeb 1995, Catelan & Theuns 1996a): 
The trend is such that smaller initial peaks, which are likely to be more asymmetric, 
are more subject to torques and have larger mean angular momenta; this would make it 
even more likely to find that low mass galaxies are also low surface brightness galaxies, 
and could potentially lead to systematic variations in galaxy-galaxy correlation properties 
as a function of surface brightness. However, given that the shift in angular momentum 
distribution with peak height has yet to be verified numerically, we will retain the simple 
assumption of the independence of angular momentum and peak height. 

Another of our assumptions is that the angular momentum of the gas is conserved 
during collapse, independent of the parameters of the collapse. This assumption should be 
robust when a collapse can be approximated as smooth; for example, in a collapse where 
there are no major mergers, or where the gas in merging subclumps is extended and hot 
enough that there is little angular momentum transfer between the gas clumps and the dark 
matter halos. However, the degree to which angular momentum conservation holds may 
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well depend on the detailed merger history of a given galaxy, which in turn could depend 
systematically on mass, angular momentum, and peak height. 

While there are theoretical and numerical tools for approximating the merger history 
of galaxy halos of a given size (e.g. Bower 1991, Bond 1991, Lacey & Cole 1993, 1994), 
there is not a consistent picture of the behavior of gas in these merging systems which 
leads to the observed galaxy structure. For example, depending on the assumptions about 
the ultraviolet background radiation and star formation feedback, galaxies in hierarchical 
formation scenarios can have from 30% to 90% of their gas cool (Navarro & White 1993, 
Evrard, Summers, & Davis 1994, Thoul & Weinberg 1996, Navarro & Steinmetz 1996). A 
consistent feature is that these simulations produce disks in which the gas is much more 
centrally concentrated than in disks of real spiral galaxies. The spin parameters of the 
simulated disks show a range of values, generally below those observed, and indicate some 
angular momentum transport to the dark matter. However, given the developing nature of 
the field and, more importantly, the failure to reproduce the observed gas profiles, it would 
imprudent to assume angular momentum transport plays a strong role in galaxy formation 
and we have no firm basis for considering its variation with the properties of galaxy halos. 

In principle, any modification can be worked into equation but, for now, we take 
the properties of the dark matter halos as representative of the properties of the resulting 
disks. This simplification may eventually prove to be a poor assumption for some individual 
galaxies, but it would be incorrect in general only if there were such extensive, violent 
merging as to entirely erase any memory of the initial angular momentum distribution from 
the resulting gas distributions, for the bulk of the galaxy population. Such a scenario seems 
unlikely as it raises serious questions about why spiral galaxies have angular momenta 
which are comparable to the predicted angular momenta of dark matter halos and how one 
avoids excessive central gas concentrations in such a model. 

In addition to the above distribution of spin angular momenta, to fully specify the 
predicted distribution of disk central surface brightnesses and scale lengths we also need 
a prediction for the distribution of galaxy masses. As a first approximation, we use a 
Schechter function with a power law of slope a at the low mass end, with an exponential 
cutoff at masses greater than M*, such that 

n(M tot ) = $*(M toi /M*) a <™ exp (-M tot /M*). (21) 

(Press & Schechter 1974, Schechter 1975). 

With the above distribution of spin angular momenta and galaxy masses, we may 
calculate the predicted distribution of disk central surface brightnesses and scale lengths. 
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First, we may use the distribution of angular momentum in equation |2(^ and the relationship 
between surface brightness and angular momentum in equation to calculate the predicted 
distribution of disk surface brightness at fixed disk mass: 



p(Z \M d ) =p(X) 



dX 



cf£ 



M d 



(22) 



which requires the following inversion of equation [TJ]: 



A(S |M d 



A + V0.944A 2 + 0.028A 
2-4A ' 



(23) 



where 
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(24) 



Then, with equations [12] and |22] and the distribution of galaxy masses n(M d ), the joint 
distribution of central surface brightness and scale length is 



n(a, S r 



n(M d , S c 



dMd 
da 



p(Z \M d )n(M d (E ,a)) 



da 



(25) 



All of the above equations involving the central surface brightness So can easily be 
reformulated in terms of magnitude per square arcsecond, fi , using the conversion 
p(^)=p(S )S /(2.51og e ). 



Because the distribution of angular momenta (eqn. |20|) is broad, spanning more than 
a factor of 10 in A, equations [H] and suggest that the distribution of galaxy surface 
brightnesses should also be broadQ. For example, keeping only the \~ 2 ~ &F dependence of 
equation |14, the log-normal distribution of angular momenta (eqn. |20D leads to a Gaussian 



Note that van der Kruit (1987) came to a nearly opposite conclusion with a similar analysis, by assuming 
that disk galaxies had one unique value of A. His Equation 15 (an analog to our equation |l^) therefore 
predicted a single characteristic disk surface density. 
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distribution of central surface brightness /i , with width <j m = 2.5(2 + 6F) log (e)cr\, or 
roughly 2.5<j\ mag/arcsec 2 . Given that most numerical studies find <3\ ~ 0.7, there is 
immediately reason to expect galactic disks to form with a very wide distribution of central 
surface brightnesses, at every luminosity. 

The peak of the (approximately) Gaussian surface brightness distribution corresponds 
to a characteristic surface brightness for each given mass: 
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(26) 



in units of B magnitudes per square arcsecond, assuming F = 0.05. The characteristic 
surface brightness for "normal" bright galaxies is therefore comparable to the canonical 
Freeman (1970) value of 21.7 ± 0.3 mag/arcsec 2 for disks. However, the expected width of 
±3.4 magnitudes suggests that there will be large numbers of bright galaxies with very faint 
disk surface brightnesses, down to ~ 24 mag/arcsec 2 ; at lower luminosities, the expected 
distribution of disk properties will shift to even lower surface brightnesses. 



3.1. The Distribution of Surface Brightness and Scale Length 

The general picture that the above sections describe is one where the structure of galaxy 



disks can be directly linked to mass and angular momentum through equations [12| and [L4[ 
At fixed angular momentum, low mass galaxies have faint central surface brightnesses and 
small scale lengths. At a fixed mass, galaxies with large angular momentum spread their 
mass over a larger area, leading to large scale lengths but faint central surface brightnesses. 
Therefore, we expect to find galaxies with a wide range of surface brightnesses and scale 
lengths, due to variations in mass and angular momenta among protogalaxies. We may 
take the distribution of surface brightness and scale length predicted by this scenario (eqn. 
5f) and compare it with the observed properties of the population of disk galaxies. 
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Fig. 4.— 

Observed and predicted distribution of exponential scale lengths a and central surface brightness 
/io (in the B band) for galactic disks. The solid circles are taken from a variety of sources 
(Romanishin et al. 1973, Boroson 1981, van der Kruit 1987, Sprayberry et al. 1995, de Blok et 
al. 1995, McGaugh & Bothun 1994, de Jong & van der Kruit 1994, de Jong 1996, Knezek 1993), 
and consist mostly of galaxies chosen from the NGC, the UGC, or the POSS-II catalog of LSBs 
(Schombert et al. 1992), all of which are either explicitly or implicitly angular diameter limited 
field surveys; however, because of how the points were assembled from the literature, the filled 
circles cannot be used to determine the relative number densities of galaxies as a function of a 
and jUQ) although they do give a good indication of the range of galaxy properties found in field 
surveys. The open squares are from the center of the Virgo cluster (Bingelli et al. 1984), the stars 
are from a catalog of LSBs in the direction of the Virgo cluster (Impey et al. 1988), including 
Malin I (upper far right). The filled triangles are the local group dwarf spheroidals compiled by 



Caldwell et al. 1992. The dotted lines are contours of constant n(lga,/xo) (from equation |25|) , 
assuming T = 3M /L0, F = 0.05, and Hq = 50km/s/Mpc (which are the same parameters 
which fit the Tully-Fisher relation in Figure ||), and assuming a\ um = —1.5, M* = 10 12 M Q , 
(A) = 0.06, and o\ = 0.7 for the mass and angular momentum distributions. The diagonal solid 
lines of positive slope are lines of constant mass (separated by factors of 10, with the heavier line 
at the characteristic mass M*), and the diagonal solid lines of negative slope are lines of constant 
angular momentum (logarithmically spaced, separated by factors of 10 ' 2 , with the heavier line at 
a characteristic angular momentum A >~ 0.06). The crossmatched region in the upper left corner 
is the area where gas pressure will prevent galaxies from collapsing (i.e. M < M m i n ; §2.2). Note 



that: 1) the region of the plane occupied by galaxies is strongly dependent on the type of survey 
used to find the galaxies; 2) cluster surveys, which are largely unbiased with regard to physical 
scale length, are dominated by galaxies with low surface brightnesses and small scale lengths, as 
would be predicted by n(lga, hq); 3) both the data and the models show the same paucity of 
physically large, high surface brightness galaxies, which reflects the exponential mass cutoff in the 
mass function of galaxies; and 4) galaxies exist with the lowest surface brightnesses to which the 
deepest surveys are sensitive, suggesting there are likely to be galaxies with even lower surface 
brightnesses than have currently been observed. See §EO for a fuller discussion. 
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In Figure ^, the distribution of surface brightness and scale length predicted by 
equation |25], n(lga,/i ), is drawn as dotted contours separated by factors of 10 in number 
density. Superimposed upon these contours are diagonal lines of constant disk mass (the 
solid lines of positive slope) and lines of constant angular momentum (the solid lines 
of negative slope). The distribution plotted in Figure [| shows several observationally 
testable features. First, the distribution is extremely broad in both surface brightness and 
scale length, suggesting that disk galaxies should form with an extremely wide range of 
properties. Second, looking at n(lg a, /io) along a line of constant mass, the distribution 
suggests that galaxies of any mass have a wide range of surface brightnesses. Thus, low 
surface brightness galaxies are not necessarily dwarf galaxies, and many examples should be 
found with masses which are comparable to normal galaxies. Third, both normal and low 
surface brightness galaxies should show the exponential cutoff with mass. At the cutoff, low 
surface brightness galaxies will have systematically larger scale lengths due to their high 
angular momentum. These effects lead to an apparently larger maximum disk scale lengths 
with decreasing surface brightness (i.e. the diaganol edge in Figure ^ towards high surface 
brightness and large scale length). Fourth, the number density rises with decreasing scale 
length and decreasing surface brightness, suggesting that there is an enormous population 
of small, low surface brightness galaxies. The rise is due to a combination of the increase 
in the galaxy mass function with decreasing mass, and of the correlation between mass and 
surface brightness and/or scale length (eqns [13] & |12|). Finally, Figure £| suggests that there 
should be a correlation between mass and surface brightness. The correlation results from 
eqn [T3], which suggests surface brightness is proportional to M 1 / 3 . However, as can be seen 
by the contours of n(\ga, fio), the correlation will be broad, due to the width of the angular 
momentum distribution. The correlation may be the source of the surface brightness - 
luminosity relation found by Bingelli et al. (1984) in the Virgo cluster, and by Ferguson & 
Sandage in Fornax (1988). 

To compare the predicted distribution of n(\ga,fio) to observations of disk galaxies, 
in Figure £| we have also plotted observed properties of galactic disks for a variety of 
surveys with different selection techniques. Because selection effects are largely responsible 
for shaping the apparent distribution, for the moment we will restrict our comparison to 
galaxies discovered exclusively in large photographic angular diamter limited field surveys 
(i.e. NGC, UGC (Nilson 1973), POSS-II LSB survey (Schombert et al. 1992), APM LSB 
survey (Impey et al. 1996)), all of which have similar angular diameter limits and thus 
similar selection criteria. These field galaxies are plotted as filled circles in Figure [|, based 
upon their published surface brightness photometry. While little can be said about the 
relative number density of the field galaxies on this plot, given that the disk properties are 
drawn from a variety sources, there are some clear trends in the range of observed disk 
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properties. First, there is an apparent cutoff at faint surface brightness, which reflects 
selection biases against LSBs; the lowest surface brightness galaxies among the filled 
circles are barely visible on the photographic plates used for discovery^. Second, there is 
an apparent cutoff at small exponential scale lengths. The cutoff results from selection 
biases in angular diameter limited surveys from which the galaxies were drawn. Angular 
diameter limited surveys are heavily biased towards finding galaxies with the largest scale 
lengths and against finding intrinsically small galaxies; unless it is incredibly close, a very 
small galaxy will not have a large enough angular extent to be included in the survey. 
Therefore, physically large galaxies will be stongly over-represented in angular diameter 
limited surveys. 

There are two other boundaries in the observed distribution of disk properties in Figure 
^ which are physically meaningful, and are not due to selection biases. First, there is an 
apparent absence of large, high surface brightness galaxies, as can be seen from the diagonal 
cutoff in the lower right. This cannot be due to selection effects, given that these galaxies 
are the easiest of all possible galaxies to detect. However, this cutoff agrees well with the 
cutoff in the predicted distribution of scale lengths and surface brightnesses, calculated 
above and plotted as the dotted contours. 

The second physically important cutoff in the distribution of disk properties is the 
absence of high surface brightness galaxies. Equation and Figure £| suggest that there 
will be significant numbers of galaxies which have surface brightnesses brighter than the 
Freeman surface brightness (1970). However, while there are many galaxies with surface 
brightnesses fainter than the Freeman value, there are no galaxies observed with surface 
brightness much brighter than the Freeman value (Allen & Shu 1979). 

In Figure ^ there is a maximum surface brightness at roughly fi Q = 19.5 mag/arcsec 2 , 
above which galaxy disks do not seem to form. We postulate that these galaxies either 
lose substantial amounts of angular momentum during collapse or become unstable to 
global modes after formation; both processes would cause the disk to loose both energy and 
angular momentum, leading to the formation of bulges and/or ellipticals. Support for this 
idea comes from a stability analysis of our model disks using the Toomre criteria for the 
stability of a rotating stellar disk (Q = 33^^ > 1 for stability (Toomre 1964), where oy r 
is the radial velocity dispersion of the stars, and k is the epicyclic frequency). While the 
Toomre criteria is derived for local disk stability, it gives a reasonable prediction of a disk's 



5 Note that historically only galaxies with (j,q < 22 mag/arcsec 2 have been well studied due to their easy 
detect ability. While this surface brightness regime spans a reasonable range in mass, it only covers a limited 
range in angular momentum, leading to the (false) impression that, not only do disks have a characteristic 
surface brightness, they have a characteristic angular momentum as well. 
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stability to the global modes which can lead to significant amounts of angular momentum 
transfer and dissipation. 

In Figure |5| we plot the Toomre parameter Q as a function of radius for galaxies with 
a range of surface brightnesses. We have assumed that all of the galaxies have the same 
central radial velocity dispersion (<jy r = lOOkm/s) and constant disk thickness, as is implied 
by observations of edge-on LSBs (Dalcanton & Shectman 1996, Cowie et al. 1996). This 
leads to a radial dependence on velocity dispersion which falls like expr/2a, as has been 
observed for normal galaxies (van der Kruit & Freeman 1986, Bottema 1993). With this 
choice of ay r , Q's inverse dependence on surface density leads to all models with surface 
densities at or below the Freeman value being stable, but leads high surface density disks, 
which are strongly self-gravitating, to be globally unstable over almost the entire disk. 

Given the large uncertainty in the velocity dispersion oy r outside of the regime of 
normal galaxies, we have also calculated the minimum velocity dispersion which would be 
needed to stabilize the disks in Figure [5[ The high surface density galaxies would require 
velocity dispersions which were typically greater than 40% of their circular velocity in order 
to be stable. However, such high velocity dispersions would lead to highly non-circular 
orbits and thus be likely to produce angular momentum transfer and dissipation. Thus, we 
again expect that the high surface density disks would not take the form predicted by our 
models. 

These very high density disks would presumably become unstable to bar formation, 
which can cause both central concentration and vertical heating, leading to the formation 
of a bulge (Combes & Sanders 1981, Pfenniger 1984, 1985, Combes et al. 1990, Pfenniger 
& Norman 1990, Friedli & Pfenniger 1990, Pfenniger & Friedli 1991, Raha et al. 1991 ). 
Although current numerical models for the secular formation of a bulge find it difficult to 
make extremely large bulges, these models usually begin with a "normal" disk, as opposed 
to an extremely high surface density disk such as the one in Figure [|; we would be most 
interested in further numerical studies of the dynamical instabilities of concentrated high 
surface density disks. Observationally, the lack of any significant color difference between 
bulges and disks, particularly in the redder bands (A(i? — K) = 0.078 ± 0.165; Peletier & 
Balcells 1996), as well as the correlation between disk scale length and bulge scale length 
(de Jong 1996, Courteau et al. 1996) and preponderance of exponential bulge profiles 
(Courteau et al. 1996) all support the notion that in a large fraction of galaxies, the creation 
of the bulge is closely tied to the formation of the disk. Of course, there are likely to be 
additional routes to bulge and/or elliptical formation; for a "messy" dissipative collapse, 
the gas may never settle into a disk at all, and entirely wind up as a condensed elliptical. 
Finally, our stability analysis shows that massive galaxies are more likely to be Toomre 
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unstable, particularly in the inner regions, which would suggest that higher mass galaxies 
should have a larger bulge-to-disk ratio, in general. Likewise, galaxies which form from low 
angular momenta halos are also more likely to be Toomre unstable, again leading to a larger 
bulge-to-disk ratio. This scenario suggests that the Hubble sequence might be a sequence 
of both mass and angular momentum. 




Fig. 5.- 

Disk stability as a function of surface brightness. The left panel plots the Toomre stability parameter 
Q = 3°3g r g S as a function of radius (scaled by the exponential scale length a) , for model disk galaxies with 
yUo = 15, 17, 19, 21, 23, &25 Smag/arcsec 2 and scale lengths of 3 kpc. The model disk with a characteristic 
Freeman value (fi ~ 21) is plotted as a dark line. Brighter disks are plotted as dashed lines, and the LSBs 
are plotted as solid lines. The models assume a radial velocity dispersion of oy T which falls off radially from 
lOOkm/s in the center such that constant disk thickness is preserved [oy r oc exp— r/2a). The high surface 
brightness disks have Q < 2 and are likely to be unstable. The LSBs, which are dark matter dominated, 
are likely to be stable. The right panel plots the minimum radial velocity dispersion needed to stabilize 
the disks in the top panel, assuming Q cr i t — 2, and plotted as a fraction of the circular velocity. The high 
surface brightness disks (dashed lines) require radial velocity dispersions which are a large fraction of their 
circular velocities in order to be stable to bar formation and other large scale modes which may result in 
angular momentum transfer and secular bulge formation. 
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The predicted distribution of scale lengths and surface brightnesses in Figure | argues 
that there should be disks over a huge range of physical scale and surface brightnesses, well 
outside of the range of disk properties represented in large catalogs such as the NGC and 
UGC. There should be large numbers of both intrinsically small galaxies and low surface 
brightness galaxies, both of which are difficult to detect. While the photographic field 
surveys which have produced the largest galaxy catalogs are basically restricted to the part 
of the a-fio plane occupied by the filled circles, other specialized surveys can find galaxies 
in other parts of the plane. For example, galaxies in the Local Group (filled triangles) and 
from deep surveys of the Virgo cluster are also plotted on Figure f|. These surveys clearly 
show that there are indeed galaxies whose properties lie well outside of the region occupied 
by galaxies from larger field surveys. More importantly, the cluster surveys plotted in 
Figure |] (open squares and stars) do not suffer from the strong bias towards intrinsically 
large galaxies which affects field surveys, and more closely reflect the intrinsic distribution 
of scale lengths; the plotted galaxies from these surveys are almost entirely found with very 
small scale lengths, suggesting that there is indeed the predicted rise towards intrinsically 
small galaxies, as shown by the contours of n(lg a, /i )- We expect the rise in the number 
of physically small galaxies to exist at both low and high surface brightnesses, suggesting 
that there may be a significant but overlooked population of very compact, high surface 
brightness galaxies; such galaxies would easily be mistaken for stars, and may well be as 
under-represented in large galaxy catalogs as the population of LSBs. 

To make the above comparison between the calculated distribution of surface density 
and the observed distribution of surface brightness, we have assumed that all disks have 
a single conversion efficiency for turning gas into stars, such that there is a single disk 
mass-to-light ratio T. In using a single mass-to-light ratio for all disks, we have swept a 
good deal of physics under the rug. In particular, we believe that our neglect of detailed 
physical processes within the galaxies causes us to systematically overestimate the surface 
brightness of LSBs, given that almost all mechanisms for reducing star formation efficiency 
are likely to be more effective in low surface density galaxies than in high surface density 
ones. First, if star formation requires a large reservoir of neutral gas, then the extragalactic 
ultraviolet background, which will completely ionize low surface density galaxies (§§), may 
suppress or shut off the formation of stars in these galaxies (Babul & Rees 1992, Efstathiou 
1992). Second, if star formation is associated with the gas in disks becoming locally Toomre 
unstable to the formation of spiral structure (implying that the Toomre stability parameter 
Q oc an/GY, is smaller than some constant, where a is the velocity dispersion of the gas, S 
is the surface density of the disk, and k is the epicyclic frequency, set by the disk structural 
parameters (see Kennicutt 1989)), then one would expect lower surface density disks to be 
less unstable for star formation. Observations by van der Hulst et al. (1993) find that LSB 
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galaxies do have HI surface densities that that fall below the critical density implied by Q 
and are about a factor of 2 lower than the HI surface densities of HSB galaxies. Third, low 
surface density, low mass galaxies are more likely to lose their gas through supernova driven 
winds than high surface density, high mass galaxies are. The gas loss both shuts off star 
formation prematurely and evolves the galaxy towards larger sizes through sudden mass 
loss and subsequent revirialization (Dekel & Silk 1986, De Young & Heckman 1994, Babul & 
Ferguson 1996) - both mechanisms which lead to lower surface brightnesses for low surface 
density objects. However, because the gas loss is expected to stop when SN-driven bubbles 
"blow-out" the top and bottom of the disk, leaving a gas void whose width is comparable 
to the thickness of the disk, these effects are likely more severe in low velocity dispersion 
halos and may not play a role in the larger LSB disks. 

In spite of these effects, for disks with central surface brightness brighter than 
/io < 24mag/arcsec 2 (for the range of masses spanned by the filled circles in Figure f|) we 
have reason to believe that the assumption of a single value of T is valid for the population 
as a whole. As discussed in § |2.1.3| , if T varied systematically with surface density over 



this range, the mean Tully-Fisher relation in Figure |3] would shift up or down towards 
brighter or fainter absolute magnitudes, as a function of surface brightness. However, 
because the mean Tully-Fisher relation is observed to be identical for all disk galaxies with 
/io < 24mag/arcsec 2 (Zwaan et al. 1995, Sprayberry et al. 1995), there cannot be significant 
variations in T with surface brightness. Likewise, because the predicted slope agrees closely 
with the measured one, there is unlikely to be a systematic variation in T with mass as 
well. Zwaan et al. (1995) do observe that the scatter in the Tully-Fisher relationship is 
about twice as high in LSBs as in normal galaxies, suggesting that individual LSB galaxies 
are more likely to have values of T which depart from the mean; however, the population 
as a whole seems to have the same mean disk mass-to-light ratio. There may well be 
systematic variations in T outside of the range of disk properties spanned by Tully-Fisher 
measurements; if included theoretically, the variation in T would be manifested in Figure f| 
as a vertical stretching of the contours of n(a,fj,o). 



3.2. Selection Effects Upon the Luminosity Function 

The predicted distribution of disk scale lengths and surface brightnesses plotted in 
Figure f| suggests that there are large numbers of uncataloged galaxies which have not 
contributed to measurements of the local luminosity function. In particular, although it is 
rarely discussed explicitly, most galaxy surveys are sensitive only to galaxies whose central 
surface brightness exceeds some minimum critical surface brightness, [i m in- As can be seen 
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from the contour lines in Figure [|, when compared to the range of surface brightnesses 
covered by major field surveys (/io < 22 mag/arcsec 2 typically), the mismeasurement of the 
luminosity function may be severe even for bright galaxies. Furthermore, because lower 
mass galaxies also tend to have lower surface brightnesses, the degree of incompleteness 
is always worse at faint luminosities; thus a surface brightness limited sample will always 
underestimate the faint end slope of the galaxy luminosity function. 




-22 -20 -18 -16 -14 

M - 5 log(h l00 ) 



Fig. 6.— 

The luminosity function which would be observed for different limiting surface brightnesses (from top to 
bottom, [iu m — oo, 24, 22, 20 mag/arcsec 2 ), appropriate for a completely unbiased survey, a deep pencil 
beam survey, a local field survey, and a multi-fiber field survey. We have assumed an intrinsic luminosity 
function of the form p(L) oc (i/i*)~ 15 exp(— L/L„), in arbitrary units of per luminosity per volume. 
Note that there is not only a significant depression of the faint end slope, but also a potentially large 
underestimate of the number of bright L* galaxies. See Figure [?] for the corresponding underestimate in the 
luminosity density. We have assumed the same parameters as in Figure 0. 
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We can quantify the degree to which the limiting surface brightness of a survey 



leads to underestimates of the luminosity function. Equation ^ gives the distribution 
of surface brightness at a fixed disk mass, M^, corresponding to a fixed disk luminosity 
L = Md/T. The observed luminosity function is therefore true luminosity function weighted 
by the integral of equation |22| from to fi m i n . Figure || shows the resulting luminosity 
function predicted by our model as a function of the surface brightness cutoff, assuming 
an underlying luminosity function n(L) oc (L/L*) -1 - 5 exp(— L/L*). The lines correspond to 
different surface brightness limits, appropriate to different types of galaxy surveys; shallow 
multi-fiber surveys of the field (e.g. LCRS (Lin 1996)), deeper surveys of the field (e.g. CfA 
(Marzke et al 1994)), or deep pencil beam redshift surveys (e.g. CFRS (Lilly et al. 1996)). 

First, from Figure |6| it is clear that the luminosity function can be significantly 
underestimated, even at L*. The field surveys with the shallowest surface brightness limits 
can be missing more than a factor of 5 of bright L* galaxies. Even deeper surveys of 
the local field can be missing a factor of 2. Only when surveys reach comparable depths 
to deep pencil beam surveys are the surveys identifying the majority of bright galaxies. 
This immediately has bearing on the apparent excess of "faint blue galaxies" at moderate 
redshifts; as pointed out by McGaugh (1994), the surface brightness selection effects are 
different in deep pencil beam surveys than in the local field galaxy surveys, potentially 
leading to the inclusion of lower surface brightness galaxies in deep surveys which are 
otherwised missed in tabulations of the local galaxy population. From Figure ^| it is clear 
that the expected magnitude of this effect is rather large, with there being 2-5 times as many 
bright galaxies visible in deep surveys as would be expected on the basis of measurements 
of the local luminosity function. We note that this is remarkably similar to the discrepancy 
between the observed luminosity density in deep pencil beam surveys and that which would 
be derived from the local field luminosity function (Dalcanton 1993). 

Second, there is a systematic decrease in the observed faint end slope with brighter 
surface brightness limits; for an intrinsic faint end slope of —1.5, the measured slope would 
be -1.15, -1.26, or -1.38 for surface brightness limits of \i m i n = 20,22, or 24. Semi-analytic 
models of galaxy formation consistently predict steeper faint end slopes than are observed, 
but without making any correction for surface brightness selection effects (Kauffmann, 
White & Guiderdoni 1993; Lacey et al. 1993; Heyl et al. 1994). While this shortfall at 
the faint end has been seen as a failure of theories of structure formation, such as CDM, 
Figure || suggests that surface brightness selection effects are easily capable of producing 
the shortfall. We should also note that while theoretical predictions of the luminosity 
function do not tend to overestimate the number of bright galaxies, this could easily be due 
to choosing normalizations based upon observations of only the population of high surface 
brightness galaxies. 
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Fig. 7, 



The differential and cumulative luminosity density in arbitrary units as a function of central surface 
brightness (upper and lower panels respectively), assuming the same parameters as in Figures || &||. Most of 
the luminosity density comes from galaxies with central surface brightnesses between 20 and 23mag/arcsec 2 
(upper panel), suggesting that current surveys are missing a significant fraction of the luminosity density 
of the universe. The lower panel shows that only deep surveys with nu m = 25 mag/arcsec 2 are capable of 
detecting most of the light in the universe. 
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a ]um =-1.5 M.= le+12 M T = 3 M /L s F = 0.05 H = 50 km/s/Mpc 
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Fig. 8.— 

The relative luminosity density as a function of B central surface brightness no and scale length a. Identical 
to Figure [|, except the dotted lines are contours of constant luminosity density, $(lga, /xo), separated by 
factors of 10. 
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The potentially large underestimate of the number density of galaxies suggests 
that current determinations of the luminosity function may also underestimate the total 
luminosity in galaxies. We can estimate the total luminosity density produced by disc 
galaxies as a function of surface brightness by integrating the intrinsic luminosity function 
weighted by p(fi \L) (equation |22|) over luminosity: 



where L min = M min /T, from equation [HJ is the expected cutoff in the galaxy luminosity 
function. Figure ^ shows the resulting luminosity density as a function of surface brightness, 
assuming a faint end slope of ai um = —1.5. The lower panel in the figure shows the 
cumulative luminosity density distribution; it suggests that only deep surveys with 
fii im ps 24mag/arcsec 2 are capable of observing most of the luminosity density of the 
universe (see, e.g., Vaisanen 1996). Figure || shows how the relative contribution to the 
luminosity density chang function of galaxy surface brightness and scale length. 

From the contours of constant luminosity density, it is clear that galaxies with central 
surface brightness as faint as 25 mag/arcsec 2 and with scale lengths as small as 0.2 kpc 
make significant contributions to the luminosity density of the universe. 

The underestimate of the local luminosity function has strong implications for 
measurements of the mass density of the universe, po- O ne common technique to measure 
Po is to take the local luminosity density and scale it by the large scale mass-to-light ratio 
measured in rich clusters such as Coma. However, if the luminosity of Coma were to be 
determined over a much wider surface brightness range than the local luminosity function, 
then the appropriate mass-to-light ratio for calculating p would be much larger than the 
measured mass-to-light ratio for Coma. Failing to correct for this effect could lead to 
measured values of po which are factors of 2 or more too low. 



We have postulated that hierarchical structure formation models naturally lead to a 
large population of low surface brightness disk galaxies. While such a pervasive population 
of LSB's is difficult to detect optically, it could potentially manifest itself in blind HI 
surveys. However, it has often been considered a failing of hierarchical structure formation 
scenarios that deep HI surveys have failed to uncover a significant population of dwarf 
galaxies. The few uncataloged dwarfs that are discovered are preferentially found near 
bright galaxies (see van Gorkom 1993 for a recent review). There does not seem to be a 
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large population of gas rich dwarfs in the small volume which has been explored in blind HI 
surveys. 

However, in light of recent work showing a sharp cutoff in HI disks at column densities 
of 10 19 cm _1 ( Corbelli, Schneider, & Salpeter 1989), the paucity of HI dwarfs may not be 
surprising. Recent work by Maloney (1993) and Corbelli & Salpeter (1993) convincingly 
demonstrates that ionization of the HI by a UV background accounts for the sharp cutoff in 
HI disks. As we have shown that low-mass galaxies tend to have low surface densities, these 
galaxies will be prone to having their hydrogen ionized, reducing their detectable HI masses 
well below their total hydrogen masses. On the other hand, HI observations of LSBs find 
that while the HI surface density of LSBs are less than normal galaxies, the ratio of HI mass 
to blue luminosity is high, suggesting that LSBs are in general less evolved than normal 
galaxies, due either to formation time, or to reduced star formation efficiency (van der Hulst 
et al. 1993, de Blok et al. 1996). This effect may partially counteract the likely decrease in 
the neutral fraction due to ionization. However, the ratio of HI mass to blue luminosity is 
only a factor of ~ 2 higher than normal galaxies, while the surface brightnesses are factors 
of ten lower, suggesting that the two competing effects are unlikely to cancel completely. 

Low surface density galaxies are therefore likely to suffer from strong biases against 
their detection in HI surveys, not just in optical surveys. A large population of LSB's could 
easily have been overlooked by existing surveys. The dwarfs that have been detected in 
HI surveys must have higher surface densities in general, and thus are more likely to have 
collapsed earlier from larger overdensities. This would explain why these dwarfs are found 
to trace the bright galaxy population. (For a more complete discussion of constraints on 
the LSB population from 21 cm surveys, see McGaugh (1996).) However, while the neutral 
gas fraction of LSBs may be small enough to keep them out of surveys of HI emission, they 
may well be seen in absorption, particularly in the Lyman-a forest. 

5. Summary 

In this paper, we have developed a scenario that links the mass and angular momentum 
of a proto-galaxy to the luminosity and surface brightness of the resultant galactic disk. 
Gas in low angular momentum proto-galaxies collapse by large factors to form high surface 
brightness galaxies. Low surface brightness galaxies form from low mass and high angular 
momentum proto-galaxies. 

As gravitational collapse in any hierarchical model with Gaussian initial conditions 
leads to a broad distribution of halo masses and angular momentum, we expect that galaxies 
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should have a wide range of surface brightness and disk scale lengths. We expect that some 
disk galaxies form with surface brightnesses less than 21 B mag/arcsec 2 , well below current 
observational limits. We also expect the number density of galaxies to rise towards small 
scale lengths and low surface brightnesses, as is seen in surveys of nearby clusters. This 
correlation also leads to an apparent correlation between surface brightness and magnitude, 
as is observed (Bingelli et al. 1984, Ferguson & Sandage 1988). 

Following Mestel (1963), Fall & Efstathiou (1980), Gunn (1982), van der Kruit (1987), 
we follow the collapse of a uniform gas cloud that has experienced a uniform external torque. 
We extend the earlier work by calculating the collapse in a realistic dark matter halo and use 
the adiabatic approximation developed by Blumenthal et al. (1986) to model the response 
of the halo to the disk. This scenario produces galaxy disks with asymptotically flat rotation 
curves and exponential light profiles, where the disk scale length is roughly proportional 
to AM 1 / 3 . Here, A is the spin parameter of the proto-galaxy and M is its mass. This 
relation leads to galaxies that obey the Tully-Fisher law as v\ oc M(r d )/r d oc M 2 / 3 oc L 2 / 3 . 
It also helps explain the disk-halo conspiracy (Bahcall & Casertano 1985) associated with 
featureless rotation curves as both the disk and halo scale lengths are proportional to M 1//3 . 

The shape of the rotation curve depends upon a galaxy's angular momentum. Low 
angular momentum disks are centrally concentrated and dominate the inner portions of the 
rotation curve. Thus, these galaxies have rapidly rising rotation curves in the centers. High 
angular momentum disks have their mass spread out to larger radii, leading to a smaller 
dynamical contribution from the disk, relative to the halo, at all radii. Thus, high angular 
momentum disks, which tend to be low surface brightness, have more slowly rising rotation 
curves. Low surface brightness disks also have apparently higher dynamical mass-to-light 
ratios, because their extended disks encompass a larger fraction of the dark matter halo 
within their optical radius. Because the disk makes a relatively small contribution to the 
galactic potential in LSBs, it serves as an excellent tracer of the initial dark mass profile and 
the shape of the dark halo. Furthermore, at a given mass, low surface brightness galaxies 
have much larger scale lengths, allowing them to probe the halo at much larger radii. 

This model predicts that very low angular momentum proto-galaxies collapse to form 
very high surface brightness galaxies. Since these very high surface brightness galaxies are 
globally unstable to non-axisymmetric perturbations, they likely form bars, bulges and/or 
ellipticals. Thus, the Toomre instability criterion may explain why few disk galaxies have 
central surface brightnesses above the canonical Freeman value. 

Our model suggests that current galaxy surveys are missing a significant fraction of 
the total number of galaxies, and have severely underestimated the faint end slope, due to 
surface brightness limits which prevent them from finding low surface brightness galaxies. 
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The luminosity density of the universe is similarly underestimated. The underestimate 
should be much less in deep pencil beam redshift surveys, and is of the right magnitude to 
explain much of the apparent "excess" in the galaxy population at moderate redshifts. It 
is intriguing that wide-field faint searches for LSBs (Dalcanton et al. 1997) are finding the 
large number of LSBs predicted by this model. 

The galaxy formation scenario makes a number of assumptions that require further 
observational and numerical testing, namely: (1) that gravitational collapse leads to a 
universal dark matter halo (Navarro et al. 1996); (2) that angular momentum is conserved 
through galaxy collapse and (3) that the initial angular momentum distribution is similar to 
that produced by a external tidal torque. All three of these assumptions are controversial. 
However, wherever possible, we have shown that our results are well supported by existing 
observations. The success of our model in matching the observations suggests that, in spite 
of the simplicity of our assumptions, the overall scenario merits careful examination. 
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7. Appendix 

In this appendix, we redo the calculation presented in the main text, but here model 
the initial halo density distribution as p(r) = Bp r /[ri(l +rf /t-q)]. This density distribution 
is close to that fit by Navarro, Frenk and White (1996) to their numerical simulations. In 
their simulations, B ~ 7 x 10 4 , and the density profile extends to r max ~ 7.4ro- The mass 
profile that corresponds to the initial dark halo mass distribution can be inverted: 



n = r 
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- 1 



(28) 



As in section 2, we can calculate the final radial distribution of the disk and the halo 
by requiring angular momentum conservation: 
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where f3 = \JG{1 - F)M tot /J max . For r max = 7.4r , the energy of the halo is 
-0.085GM^/r , thus, r = 0.141a, where a is the radius of the initial fluctuation that 
collapsed to form the galaxy. Hence, (3 ~ 0.15/A. Figure |9| shows the rotation curves and 
surface density profiles for a range of values for A and F. For this model, the disk scale 
length, r disk ~ a\ b r , where a = 4.6 * (1 + 20/) and 6=1 + 6/ 

Following the same approach used in sections 2 and 3, we can compute the predicted 
distribution of surface densities at a given luminosity: 
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This model yields similar predictions to the Hernquist model used in the main text; 
however, now the width of the central surface density distribution is slightly broader: 
^ = 1.08*5(1 + 6/). 




Fig. 9.— 

Same as Figure [l], but for an initial halo density profile of the form: p(rj) oc r^(rf + Tq)~ 



